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Introduction 


Ceramic  materials  are  increasingly  being  used  as  armor  elements,  engine  turbine 
blades,  and  other  structural  elements  because  of  their  enhanced  dynamic  compressive 
strength  and  high  temperature  properties.  Understanding  of  ceramic  behavior  under 
impact  loading  conditions  is  essential  in  the  design  of  improved  impact  resistant  materials 
for  dynamic  structural  and  armor  applications.  For  this  purpose,  a  few  impact/shock 
wave  propagation  experimental  techniques  are  available.  The  experimental  data  are  often 
used  in  the  development  of  constitutive/damage  models  to  describe  the  complex  impact 
behavior  of  ceramic  materials.  An  increasing  use  of  advanced  Finite  element/difference 
computer  codes  in  the  armor/anti-armor  design  analyses  requires  development  of 
dependable  and  realistic  material  models  in  the  codes.  This  requirement  is  all  the  more 
important  for  advanced  materials  such  as  the  ceramics  and  glass  Fiber-reinforced  plastic 
materials. 

Recently,  Rajendran  (1,2)  presented  internal  state  variable-based  constitutive  equations 
to  describe  the  impact  behavior  of  ceramic  materials.  The  model  constants  for  AD85 
ceramic  were  determined  using  the  data  from  a  combination  of  impact  experiments  includ¬ 
ing  the  uniaxial  strain-based  plate  impact  experiment  and  the  uniaxial  stress-based  split 
Hopkinson  bar  (SHB),  as  well  as  the  rod-on-rod  experiments  (3,4).  The  ceramic  constitu¬ 
tive  model  assumes  the  following: 

•  Pre-existing  randomly  distributed  flaws  in  the  ceramics 

•  Plastic  flow  in  the  ceramics  when  shocked  about  the  HEL 

•  No  plastic  flow  in  tension 

•  Degradation  of  elastic  moduli  under  both  compression  and  tension 

•  Pulverization  occurs  under  compressive  loading  when  the  accumulated 
crack  density  reaches  a  value  of  0.75.  The  flaws  are  activated  according 
to  a  generalized  GrifFith  criterion 

There  are  three  constants  that  characterize  the  virgin  ceramics;  number  of  micro- 
flaws  (N*),  initial  size  of  the  flaw  (ao),  and  dynamic  coefficient  of  friction  (|i).  For 
damage  evolution,  the  model  has  two  constants;  crack  growth  factor  (ni)  and  crack 
growth  index  (n2).  These  two  loading-state  dependent  constants  take  different  values 
under  tension  and  compression.  For  most  ceramics,  ni  and  nj  were  assumed  to  be  equal 
to  one  under  tensile  crack  growth.  The  mode  I  fracture  toughness  (Kic)  is  determined 
through  direct  measurements  or  from  handbooks.  At  present,  representative  values  for 
the  model  parameters  were  determined  based  on  the  model’s  ability  to  reproduce  the 
measured  stress/velocity  proFiles. 

This  work  describes  the  modeling  of  three  armor  ceramics:  titanium  diboride 
(TiB2),  silicon  carbide  (SiC),  and  boron  carbide  (B4C).  The  experimental  data  of  Kipp 
and  Grady  (5)  were  employed  in  the  modeling  efforts.  Steinberg  (6),  Johnson  and 
Holmquist  (7),  and  Addessio  and  Johnson  (8)  also  employed  the  Kipp  and  Grady  data 
in  their  modeling  efforts.  Recently,  Mandell  and  Henninger  (9)  evaluated  the  ceramic 
models  of  Steinberg  and  Johnson-Holmquist  using  the  MESA  code  (10). 
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In  the  Description  of  the  Model  Section  the  salient  features  of  the  ceramic  model 
given  in  References  1  and  2  are  briefly  described.  The  Results  and  Discussion  Section 
shows  capabilities  of  the  model  in  reproducing  the  VISAR  data  from  the  plate  impact 
experiments.  The  various  features  of  the  ceramic  model  are  explored  through  a  series  of 
simulations  in  which  the  model  constants  were  systematically  varied.  The  damaged 
moduli  expressions  of  Margolin  (11)  are  provided  in  the  Appendix. 


Description  of  the  Modei 

The  ceramic  model  assumes  the  existence  of  flaws  (microcracks)  in  the  virgin 
material.  Therefore,  the  model  does  not  require  any  microcrack  nucleation  criterion.  The 
existing  flaws  are  activated  (extended)  when  the  stress  state  satisfies  the  generalized 
Griffith’s  criterion.  The  effect  of  crack  interaction  due  to  coalescence  is  neglected. 
When  the  accumulating  crack  density  reached  a  critical  value,  the  model  assumes 
pulverization  of  the  ceramic. 

Constitutive  Relationships 

The  ceramic  model  assumes  strains  due  to; 


•  Elastic  and  dislocation-based  plastic  deformations 

•  Microcrack  opening  and  sliding 

•  Pore  collapsing 


In  this  section,  the  stress-strain  relationship  for  a  microcracked  material  is  briefly 
discussed.  The  total  strain  is  decomposed  into  elastic  ( £  )  and  plastic  strains  (  £  )  as 


+ 


(1) 


where  the  elastic  strain  consists  of  the  elastic  strain  of  the  intact  matrix  material  and 
the  strain  due  to  crack  opening/sliding.  The  plastic  strains  associated  with  pore  collaps¬ 
ing  are  modeled  using  the  equations  derived  by  Rajendran,  et  al.,  (12)  for  porous  ductile 
aggregate.  Plastic  flow  and  pore  collapse  are  assumed  to  occur  in  the  ceramic  only 
when  the  applied  pressure  exceeds  the  pressure  at  the  Hugoniot  elastic  limit. 

The  elastic  stress-strain  equations  for  the  microcracked  aggregate  material  are  given  by: 


(2) 


The  components  of  the  stiffness  tensor  M  following  Margolin’s  work  (11)  are  given  in 
Appendix  A.  The  total  stress  Ojj  is  decomposed  into  deviatoric  stress  (Sy)  components 
and  pressure  (P): 


(3) 


The  pressure  is  calculated  through  the  Mie-Gruniesen  equation-of-state  which  is  given  by: 


p  =  [PH(i-o.5rp)  +  rp^(i-ijj  (4) 

where 

Ph  =  Ky  (  Pi  Ti  +  |32  Ti^  +  p3  )  .  (5) 

Ky  is  the  ratio  of  microcracks-degraded  bulk  modulus  and  intact  bulk  modulus. 

This  ratio  reduces  from  1  to  0  as  the  elastic  properties  degrade  with  microcracking. 

T)  (=  -  I  )  is  the  elastic  volume  strain  (ey)  and  P|,  P2,  and  P3  are  the  empirical 

parameters.  G  is  the  Mie-Gruneisen  parameter,  Po  is  the  material’s  initial  density,  lo  is 
the  initial  value  of  internal  energy,  and  I  is  the  current  internal  energy. 

The  elastic  deviatoric  stress-strain  relationship  is  given  by: 

S..  =  2G  (e..  -  e.P)  (6) 

ij  ij  ij 

where  ejj  are  the  total  deviatoric  strains  and  G  is  the  degraded  shear  modulus.  The 
intact  bulk  and  shear  moduli  were  corrected  for  porosity  using  the  correction  factors  as 
suggested  by  Johnson  (13).  The  initial  porosity,  fo,  is  a  material  model  constant.  The 
porosity  is  assumed  to  decrease  due  to  pore  collapsing  at  pressures  above  the  HEL. 

Wh'*.n  the  estimated  von-Mises  stress  ( V  3  J^  )  exceeded  the  compressive  strength  (Y) 
measured  from  the  plate  impact  experiments,  plastic  flow  is  assumed  and  the  deviatoric 

stresses  were  determined  by  satisfying  the  von-Mises  criterion  ( Y  =  ^  3  )■  The  strain 

rate  dependent  compressive  strength  was  described  by  a  simple  relationship: 


Y  =  Cl  (1  -I-  Csin  e) 


(7) 


where  Ci  and  C3  are  model  constants  and  e  is  the  applied  strain  rate.  The  constant, 
Cl,  is  assumed  to  be  equal  to  Yhel  which  is  calculated  using  the  formula: 


C,  =  Y 


HEL 


HEL 


(—  +1) 
2G 


(8) 


Ohel  is  the  maximum  amplitude  of  the  elastic  shock.  The  constant,  C3,  is  calibrated  by 
matching  the  "plastic  wave"  portion  of  the  stress  or  velocity  profiles  of  conventional 
planar  plate  impact  experiments. 
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Definition  of  Damage 

In  the  ceramic  model,  microcrack  damage  is  measured  in  terms  of  the  dimensionless 
microcrack  density  y  which  is  defined  as; 

Y  =  N*  a^  •  (9) 

N*  is  the  average  number  of  microflaws  per  unit  volume  and  a  is  the  maximum  mi¬ 
crocrack  size  which  is  treated  as  an  internal  state  variable.  The  initial  values  of  these 
two  parameters  are  material  model  constants.  The  microcrack  is  assumed  to  extend  when 
the  stress  state  satisfies  the  Griffith’s  criterion  (14,15).  The  extension  of  microcracks 
causes  stress  relaxation. 


The  Griffith  criterion  (14),  employed  in  the  model,  includes  the  fracture  toughness 
Kic,  as  well  as  a  dynamic  frictional  coefficient  |i  as  model  parameters.  Since  N*  is 
assumed  to  be  a  constant  in  the  model,  the  increase  in  y  entirely  due  to  increase  m 
the  crack  size.  Experimental  measurements  to  determine  the  initial  number  of  microflaws 
in  a  unit  volume  and  the  initial  maximum  size  of  the  microcracks  are  very  difficult  and 
tedious.  The  values  are  usually  guessed  and  calibrated  based  on  the  model’s  ability 
to  reproduce  the  measured  stress/velocity  histories  in  plate  impact  experiments. 

Damage  Growth  Model 


The  crack  extension  (damage  evolution)  law  is  derived  from  a  fracture  mechanics- 
based  relationship  (16)  for  a  single  crack  propagation  under  dynamic  loading  conditions 
and  the  law  is  described  by: 


r 

a  =  ni  Cr 


"2  -] 

'  Gc^l  I  Gi<Gc,a=0 
G|  J  (ji  >  Gc  ,  a  >  0 


(10) 


where  Cr  is  the  Rayleigh  wave  speed,  Gc  is  the  critical  strain  energy  release  rate  for 
microcrack  growth,  and  Gi  is  the  applied  strain  energy  release  rate.  The  n|  is  used  to 
limit  the  microcrack  growth  rate.  The  crack  growth  index  n2  affects  the  crack  growth 
Gc 

rate  when  the  ratio  -  is  less  than  1,  but  close  to  1.  Different  crack  extension  rates 

G, 

are  assumed  for  tensile  and  compressive  loading  conditions  and,  therefore,  n]  and  n2 
are  different  (ni  is  n“  and  n2  is  n“  for  crack  sliding;  ni  is  n|  and  n2  is  n^  for 
crack  opening,  accordingly).  The  ceramic  is  assumed  to  pulverize  under  compression 
when  y  reached  a  critical  value  of  0.75  as  described  in  Reference  1.  The  strength  of  the 
post-fractured  ceramic  is  described  by: 


Y  =  a  +  P  P  (II) 

where  a  and  P  are  the  model  constants.  In  the  present  work,  values  for  these  two  post¬ 
fracture  model  constants  are  arbitrarily  assumed  to  be  equal  to  0  and  1,  respectively.  In 
a  later  section,  the  effect  of  P  on*  the  stress  computed  profiles  are  presented  to  justify 
the  assumption  of  one  for  p. 
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Ptate  Impact  Experiment 

Plate  impact  test.>  provide  a  loading  path  that  is  very  different  from  conventional 
SHB  tests.  The  '"^formation  isr  that  of  one-dimensional  strain,  and  the  mean  stress  is 
generally  very  high  compared  to  that  in  the  SHB  tests  under  one-dimensional  stress  state. 
Strain  rates  are  10^/s  or  higher.  The  material  undergoes  compression  followed  by  ten¬ 
sion.  In  general,  plate  impact  experiments  are  essential  for  calibrating  and  validating 
high  strain  rate  material  models  that  aspire  to  general  applicability. 

Figure  1  shows  typical  velocity  profiles  that  correspond  to  different  levels  of  spalla¬ 
tion,  including  a  case  where  spallation  does  not  occur.  The  challenging  part  of  modeling 
the  impact  behavior  of  ceramics  is  in  the  model’s  ability  to  reproduce  these  profiles  at 
various  velocity  levels  using  the  same  set  of  model  parameters.  The  various  levels  of 
spall  occur  according  to  the  impact  velocity  levels  for  a  given  flyer-target  configuration. 
The  thickness  and  material  of  the  flyer  plate  affects  the  pulse  duration.  Model  constants 
were  determined  by  matching  the  computed  velocity  profiles  with  experimentally  measured 
profiles.  By  trial  and  error,  the  model  constants  are  improved  from  the  initial  guess. 

A  set  of  constants  based  on  the  best  match  between  the  simulated  time  histories  of  stress 
and/or  velocity  and  the  experimental  data  is  the  preliminary  values  for  the  constants. 

This  set  may  not  assure  the  generality  of  the  model  constants;  however,  a  few  constants 
can  be  further  improved  by  calibrating  them  against  data  from  some  other  experimental 
configurations.  Rajendran  and  Grove  (4)  exercised  this  approach  for  AD8S  ceramic  and 
determined  a  set  of  constants  that  were  configuration  independent. 
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The  plate  impact  tests  considered  in  the  present  study  were  conducted  by  Kipp  and 
Grady  (5)  to  determine; 


•  Hugoniot  elastic  limit 

•  Spall  threshold  of  SiC,  B4C,  and  TiB2 

•  Stress-strain  path. 

Kipp  and  Grady  provide  a  detailed  discussion  on  their  plate  impact  experiments  in 
Reference  S.  The  model  parameters  determination  scheme  fully  utilized  the  plate  impact 
test  data  of  Kipp  and  Grady  for  SiC,  B4C,  and  TiB2. 


Results  and  Analyses 

This  section  describes  the  shock  response  to  the  three  ceramic  materials  based  on  the 
plate  impact  experimental  simulations  using  the  ceramic  model  presented  in  the  previous 
section.  In  the  simulation,  the  ceramic  strength  and  stiffness  were  assumed  to  degrade 
due  to  microcracking.  Recently,  Grove  (17)  implemented  the  ceramic  model  into  the 
1991  version  of  the  EPIC  code  (18).  To  determine  the  constants  for  SiC,  B4C,  and 
TiB2  using  the  VISAR  data  of  Kipp  and  Grady  (5),  this  version  of  the  EPIC  code  was 
employed.  The  dimensions  of  the  flyer  and  target  plates,  as  well  as  the  impact  velocities 
in  the  plate  impact  experiments  of  Kipp  and  Grady,  are  given  in  Table  1.  The  shock 
stresses  (maximum  axial  stress)  in  these  experiments  were  two  to  four  times  the  Hugoniot 
elastic  limit  of  the  ceramics.  The  window  material  was  lithium  fluoride  whose  density 
was  2641  kg/m^. 


Table  1 .  Kipp  and  Grady  (5)  plate  impact  expetimental  details 


Test  no. 

Material 

Impact  velocity 
(m/s) 

Ryer  thickness 
(mm) 

Target  thickness 
(mm) 

1 

SiC 

1542 

3.987 

8.939 

2 

SiC 

2100 

3.995 

8.940 

3 

B4C 

1546 

3.920 

9.044 

4 

B4C 

2210 

3.917 

9.033 

5 

TiBa 

1515 

3.972 

10.804 

6 

TIB2 

2113 

3.337 

10.747 

Kipp  and  Grady  performed  plate  impact  experiments  at  two  velocities  for  each  of  the 
three  ceramics;  one  at  a  lower  level  of  about  1.5  km/s,  and  the  other  at  a  higher  level 
of  about  2.1  to  2.2  km/s.  The  ceramic  model  constants  were  determined  using  the  lower 
velocity  experiments.  In  other  words;  only  tests  ',  3,  and  5  were  employed  in  the  model 
constants  calibration  scheme.  Later,  the  VISAR  results  of  the  higher  velocity  tests  2,  4, 
and  6  were  predicted  through  EPIC  simulations  by  employing  the  values  determined  from 
the  low  velocity  tests. 
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Mesh  and  Time  Step  Sensitivity  Study 

Since  the  governing  constitutive  equations  are  coupled,  nonlinear  equations,  a  stable 
and  accurate  solution  demands  a  robust  numerical  scheme  whith  will  assure  not  only  con¬ 
vergence  of  the  solution  but  also  yield  mesh  and  time  step  independent  solutions.  Grove 
(17)  used  a  diagonally  implicit  Runge-Kutta  scheme  in  the  implementation  of  the  model 
into  the  EPIC  code.  The  simulations  utilized  the  one-dimensional  strain  option  in  the 
EPIC  code  to  model  the  plate  impact  configuration.  To  investigate  the  effect  of  time 
step  and  mesh  size  on  the  numerical  results,  two  different  time  steps  and  meshes 
were  respectively  considered.  A  base  line  profile  for  a  0.2  mm  mesh  was  obtained 
from  a  plate  impact  simulation.  The  model  parameters  used  in  the  base  line  profile  are 
N*  =  10^; a©  =  0.005mm;  Kic  =  4 MPa  Vnr ;  |J.  =  0.1;  n~  =  0.1  .  In  the  sensitivity  study, 
the  effect  of  mesh,  time  step,  and  the  model  parameters  on  the  simulated  velocity  profile 
is  compared  with  this  base  line  profile. 

To  investigate  the  mesh-size  effect,  simulation  with  a  0.1  mm  mesh  (fine  mesh)  was 
performed  and  the  corresponding  velocity  profile  is  compared  with  the  base  line  profile 
(solid  line)  in  Figure  2.  The  results  were  almost  identical,  as  can  be  seen  from  Figure 
2.  The  kink  at  the  HEL  in  fine  mesh  was  more  pronounced  than  the  corresponding  kink 
in  the  coarse  mesh.  The  time  step  effects  were  studied  using  two  time  steps;  one 
with  a  maximum  allowable  time  step  of  0.01  ms  and  two  with  0.001  ms.  The  velocity 
time  histories  from  these  two  simulations  were  almost  identical,  as  can  be  seen  from 
Figure  3,  confirming  minimum  influence  of  the  time  step  on  the  numerical  results. 


Figure  2.  The  effect  of  mesh  size  on  the  computed  velocity  profile. 


7 


TIME  STEP  EFFECT 


Figure  3.  The  effect  of  time  step  on  the  computed  velocity  profile. 


Sensitivity  Study  on  Modei  Constants 

There  are  eight  constants  in  the  ceramic  model  to  describe  the  microcracking 
behavior:  N*,  a©,  M-,  n~  n“  n|,  n^,  Kjc  .  For  tensile  loading  (negative  pressure),  the 
crack  growth  factor  n|  and  the  crack  growth  index  n^  are  assumed  to  be  equal  to 
one.  Under  compressive  loading  (positive  pressure),  the  best  suitable  values  for  n~ 
and  n“  are  determined  by  matching  the  computed  profile  with  the  plate  impact  experi¬ 
mental  data.  To  understand  the  significance  of  several  of  the  model  constants  on  the 
shape  and  amplitude  of  the  measured  velocity  histories,  a  sensitivity  study  was  carried 
out  in  the  present  work.  Direct  measurements  of  the  microcrack  parameters  ao,  N  * , 
and  p.  (dynamic  coefficient  of  friction)  are  extremely  difficult.  The  estimation  of 
realistic  values  for  these  constants  requires  microstructural  characterization  of  the 
virgin  ceramic  through  both  nondestructive  and  destructive  test  techniques. 

In  the  present  modeling  efforts,  due  to  a  lack  of  experimental  measurements,  the 
initial  microcrack  size  is  assumed  to  be  equal  to  a  fraction  of  average  grain  size;  there¬ 
fore,  the  number  of  flaws  is  determined  through  numerical  simulations  of  the  plate  impact 
experiments.  In  effect,  there  are  only  four  adjustable  constants  to  describe  the  microcracking 
behavior:  N  * ,  ao,  p.,  and  n  ~ .  Therefore,  the  main  idea  behind  the  sensitivity  study  is 
to  identify  a  range  of  values  for  these  four  model  parameters  and^  isolate  the  effect  of 
any  one  parameter,  or  a  combination  of  parameters,  on  certain  portions  of  the  stress/velocity 
history.  Based  on  the  sensitivity  analyses,  it  may  then  be  possible  to  develop  a  stand¬ 
ard  methodology  to  estimate  appropriate  values  for  the  model  constants. 
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Number  of  Flaws,  N  * 

Three  simulations  were  considered  using  the  values  of  lo",  10^  (base  line  case),  and 
10^  for  N*.  The  corresponding  results  are  presented  in  Figure  4.  A  value  of  10^ 
did  not  produce  complete  spall  fracture.  The  ceramic  retained  some  of  its  spall  strength, 
as  one  can  see  from  the  dashed  line  curve  beyond  point  "S".  The  intermediate  value  10^ 
caused  partial  compressive  cracking  followed  by  tensile  cracking.  When  the  compressive 
microcracking  is  absent  for  N*  =  10^,  the  pulse  duration  is  according  to  the  intact 
ceramic  properties.  In  this  case,  the  microcracking  of  the  ceramic  begins  during 
unloading  and  spalls  the  ceramic  completely,  as  indicated  by  the  profile  between  points 
"D"  and  "E"  of  the  solid  curve.  Based  on  these  simulations,  a  value  of  10^  seems  to  be 
a  possible  value  for  in  order  to  cause  spall  without  excessive  compressive  damage. 
A  high  value  of  10**  caused  excessive  damage  under  compression  and  complete  spall 
occurred.  The  velocity  profile  between  point  "A"  and  point  "C"  significantly  differs  from 
the  results  for  10^  and  10^;  this  can  be  explained  by  analyzing  the  pressure,  damage, 
and  axial  stress  histories  inside  the  target. 

Effect  of  Number  of  Flowe 


Figure  4.  The  effect  of  the  number  of  flaws  on  the  computed 
velocity  profile. 


For  this  purpose,  an  element  very  close  to  the  impact  plane  is  considered  and  the 
time  histories  of  these  three  quantities  are  plotted  in  Figure  S.  The  pressure  and  stress 
are  normalized  with  respect  to  ohel  so  that  the  damage  history  can  be  presented  in  the 
same  plot.  A  value  of  15  GPa  was  used  in  the  normalization.  The  microcracking 
behind  the  shock  wave  under  compressive  pressure  (the  code  treats  the  compressive  pres¬ 
sure  as  a  positive  quantity)  causes  damage  in  the  ceramic.  The  damage  reaches  one  at 
about  0.2S  |i,s. 
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PRESSURE.  STRESS.  AM)  DAMAGE  HISTORIC 
NEAR  THE  irf>ACT  PUtfC  INSIDE  T)£  TARGET 


Figure  5.  The  time  histories  of  pressure,  stress,  and  damage' 
in  an  element  near  the  impact  plane 


Upon  high  velocity  impact  loading  conditions,  the  ceramic  material  pulverizes  under 
compressive  pressure.  In  Figure  5,  pulverization  occurs  at  about  3  p.s.  The  normalized 
pressure  level  continues  to  remain  at  the  same  level  beyond  the  point  of  pulverization. 
Since  the  bulk  modulus  does  not  degrade  under  compressive  pressure  according  to  the 
model,  especially  when  the  three  principal  stresses  are  compressive,  the  shock  pressure 
level  did  not  change. 

Due  to  absence  of  dilatation  under  high  confinement  pressure,  the  density  of  the 
confined  granular  ceramic  remains  the  same  as  that  of  the  intact  ceramic.  When  one  of 
the  three  principal  stresses  becomes  tensile  while  the  pressure  remains  under  compression, 
crack  opening  could  occur  causing  dilatation.  Under  such  loading  conditions  the  bulk 
modulus  will  degrade  in  the  modeling.  This  is  not  the  case  for  an  element  that  is  close 
to  the  impact  plane,  as  can  be  seen  from  the  normalized  (principal)  stress  plot  (see  the 
dashed  line  in  Figure  5).  This  axial  stress  remains  compressive  eliminating  any  bulk 
modulus  degradation  possibilities. 

Initial  Microcrack  Size,  ao 

To  investigate  the  effect  of  microcrack  size  on  ceramic  fracture,  three  simulations 
with  values  0.1  micron  (|i.m),  0.2  pm,  and  0.5  |im  for  a©  were  considered.  In  these 
simulations,  p.  was  0.1  and  N*  was  10^.  The  velocity  time  histories  for  the  three 
values  are  plotted  in  Figure  6.  A  small  flaw  size  of  0.1  pm  and  below  did  not  cause 
any  microcracking. 
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Initial  Crock  Length  Effect 


Figure  6.  The  effect  of  initiai  crack  size  on  the  computed 
velocity  profile. 


An  additional  simulation  with  ao  =  0  reproduced  the  result  of  the  0.1  pm  crack. 

The  0.2  pm  flaw  caused  partial  spall.  The  crack  sizes  greater  than  0.3S  pm  induced 
complete  spall  in  the  target.  A  microcrack  size  between  the  range  of  0  to  0.35  did  not 
influence  the  loading  wave  portion  of  the  profile  between  points  "A"  and  "C"  as  well  as 
the  elastic  release  between  points  "C"  and  "S".  The  range  will  vary  according  to  differ¬ 
ent  values  of  the  other  two  microcrack  parameters,  p  and  N  * . 

Coefficient  of  Friction,  p 

Figure  7  shows  the  effects  of  p  on  the  velocity  profile  for  n  “  =  0.2 .  The  main  in¬ 
fluence  of  this  parameter  is  to  decrease  the  effective  shear  stress  (x  -  pP)  level  to  sat¬ 
isfy  the  Giffith  criterion  under  compressive  pressure,  x  is  the  applied  shear  stress  and  P 
is  the  pressure  on  the  crack  surface.  Several  simulations  were  performed  to  investigate 
this  influence.  The  results  indicated  that  the  role  of  this  parameter  was  closely  tied  to 
the  values  chosen  for  the  other  microcrack  constants,  especially  the  compressive  crack 
growth  factor  n~.  For  instance,  the  effect  of  p  was  insignificant  for  n~  =0.1; 
whereas,  the  effect  became  significant  for  n“  =  0.2 
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EFFECT  OF  COEFFICIENT  OF  FRICTION 


Figure  7.  The  effect  of  coefficient  of  friction  on  the  computer 
velocity  profile. 


The  lower  values  caused  premature  microcracking  while  the  shock  amplitude  was 
still  increasing  during  the  "plastic"  wave  ramp  between  points  "B"  and  "C".  Since  the 
compressive  microcracking  occurred  well  above  the  threshold  level  of  "complete  spall 
(under  tension),"  the  influence  of  on  the  spall  signal  was  negligible.  This  study  sug¬ 
gests  that  a  minimum  value  which  will  not  cause  any  premature  compressive  damage 
shall  be  assigned  to  p..  In  fact,  this  minimum  value  will  allow  the  independent  deter¬ 
mination  of  other  constants.  In  the  modeling,  a  value  of  0.1  was  assumed  for  p  for 
SiC,  B4C,  and  TiB2. 

Compressive  Crack  Growth  Factor,  n  ~ 

This  constant  limits  the  damage  growth  rate  (see  Equation  10)  under  compression  to 
njCR.  When  the  Griffith  criterion  is  satisfied  for  a  given  stress  state,  the  time  dependent 
microcrack  extension  controlled  the  rate  at  which  the  elastic  moduli  are  degraded.  This 
degradation  affects  the  shock  wave  speed  and  this,  in  turn,  affects  the  pulse  duration. 
Experimentally  measured  values  are  available  for  the  limiting  dynamic  crack  propaga¬ 
tion  velocity  under  the  crack  opening  mode.  However,  such  measurements  are  not  avail¬ 
able  for  dynamic  crack  propagation  under  mode  II  (crack  sliding  under  shear).  The  n  ~ 
effect  on  the  computed  velocity  profile  is  shown  in  Figure  8. 

A  value  that  is  lower  than  0.1  did  not  introduce  any  further  change  in  the  profile. 
The  increasing  values  between  0.1  and  0.2  widened  the  pulse  width  due  to  elastic  moduli 
degradation.  TJie  shock  amplitude  diminished  as  the  microcracking  increased  with 
increasing  ni.  The  spall  rebound  portion  of  the  velocity  profile  was  not  significantly 
affected  when  n  ~  was  varied  in  the  simulations.  This  is  because  the  degradation  of 
the  elastic  properties  during  compression  is  well  above  the  threshold  to  cause  com¬ 
plete  spall  upon  the  arrival  of  tensile  stress  wave. 
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EFFECT  OF  COFFRESSIVE  CRACK  GROUTH  FACTOR 


TIME  (MICROSECONDS) 


Figure  8.  The  effect  of  compressive  crack  growth  factor  on  the 
calculated  velocity  profile. 


Compressive  Crack  Growth-  Index,  n  ^ 

To  verify  the  effect  of  n“  on  the  velocity  profile,  several  simulations  were  per¬ 
formed.  The  results  showed  that  the  computed  profile  was  unaffected  by  n  “ .  Since 
the  ratio  between  Gc  and  Gi  is  large  at  impact  velocities  above  1  km/s,  the  values  of 
this  index  between  0.001  and  1  did  not  alter  the  ceramic  response.  However,  two  addi¬ 
tional  simulations  at  a  lower  velocity  also  produced  very  similar  wave  profiles.  In  these 
simulations  the  initial  microcrack  size  was  assumed  to  be  equal  to  I  |am.  Grove,  et  al.,  (3) 
found  that  this  index  influenced  the  computed  stress  histories  in  a  ceramic  rod  due  to  an 
impact  by  a  similar,  but  shorter  rod  at  low  velocity.  In  their  study,  the  states  of 
stress  and  strain  at  the  gauge  location  were  one-dimensional  (uniaxial)  stress  and  three-di¬ 
mensional  strain,  respectively,  whereas  in  the  plate  impact  configuration,  the  stress  state 
was  three-dimensional  and  the  strain  state  was  one-dimensional. 

Tensile  Crack  Growth  Factor  n  | 

Ceramics  are  weak  under  tensile  loading.  The  fracture  toughness  of  a  ceramic  is 
very  low  compared  to  metals.  Upon  high  velocity  impact,  the  compressive  shock  loading 
causes  damage  in  the  ceramics.  When  the  tensile  wave  arrives  at  locations  where  the 
ceramic  has  been  partially  damaged  due  to  the  initial  compressive  shock,  it  induces  ad¬ 
ditional  damage  leading  to  complete  spall.  The  computed  velocity  profiles  for  n  |  equal 
to  0.2,  0.3,  and  1.0  are  shown  in  Figure  9. 
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EFFECT  OF  TENSILE  CRACK  GROWTH  FACTOR 


Figure9.  The  effect  of  tensle  crack  growth  factor  on  the 
computed  velocity  prof9e. 


The  crack  density  levels  increase  with  an  increasing  value  of  n|.  Since  the  degrada¬ 
tion  of  spall  strength  is  not  complete  for  lower  values  of  n|,  the  calculated  profiles 
showed  pullback  signals  (spall  rebounds)  indicating  nonzero  spall  strength.  The  value 
of  1  caused  complete  spall  in  the  ceramic  target,  as  can  be  seen  from  the  velocity  pro¬ 
file  between  points  "E"  and  "F"  (solid  curve).  By  comparing  these  simulations  with 
an  experimental  profile,  a  suitable  value  for  tensile  ni  can  be  chosen  while  modeling 
the  impact  behavior  of  any  ceramic.  Experimental  measurements  on  a  single  crack 
propagation  under  dynamic  loading  conditions  suggest  that  the  Rayleigh’s  wave  speed 
is  the  limiting  value  for  this  factor;  therefore,  in  the  present  analysis  a  value  of  I  was 
used  for  all  three  ceramics  eliminating  the  need  for  determining  the  constants  from  the 
plate  impact  experiments. 

Tensile  Crack  Growth  Index,  n  ^ 

This  parameter  appears  in  the  damage  evolution  law  of  Equation  10  as  02.  The 
ratio  of  critical  strain  energy  release  rate  (with  respect  to  crack  extension),  Gc,  and 
applied  strain  energy  release  rate,  G|,  is  modified  through  n2.  The  ratio  is  equal  to  1 
when  Gc  =  Gi.  The  crack  growth  occurs  when  the  ratio  is  below  1.  As  soon  as  the 
Griffith  criterion  is  satisfied  and  when  the  ratio  is  closer  to  1,  this  index  influences  the 
microcracking  rate.  However,  the  effect  on  the  computed  velocity  profile  is  significant 
only  for  extremely  low  values  of  the  index  as  shown  in  Figure  10. 
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EFFECT  OF  TENSILE  CRACK  C3ROUTH  INCEX 


Figure  10.  The  effect  of  tensile  crack  growth  index  on  the 
computed  vekx:ity  profile. 

A  value  of  0.01  suppressed  the  spall  fracture,  though  the  Griffith  criterion  was  satis¬ 
fied.  The  larger  values  0.1  and  1.0  produced  complete  spall  and  the  corresponding  pro¬ 
files  were  similar;  therefore,  to  cause  complete  spall,  a  value  between  0.1  and  I  can  be 
chosen  for  the  tensile  crack  growth  index  n2.  For  simplicity,  a  value  of  1  is  chosen  for 
the  three  ceramics. 

Fracture  Toughness,  Kic 

The  ceramic  model  does  not  contain  any  microflaws  nucleation  model.  A  number  of 
flaws  are  assumed  to  exist  in  the  virgin  ceramic  materials.  These  flaws  are  activated 
when  the  stress  state  satisfies  the  Griffith  criterion  for  an  initial  maximum  flaw  of  size, 
ao.  The  value  of  Kic  between  1  MPaVjjT  and  20  MPaVitT  did  not  show  any  effect  on  the 
velocity  history.  In  the  modeling,  the  static  fracture  toughness  values  from  handbooks 
were  used  for  the  three  ceramics. 

Post  Fracture  Strength  Parameter  p 

The  plate  impact  experimental  configuration  seems  to  be  insensitive  to  this  parameter. 
The  simulations  with  three  different  values  for  P  showed  hardly  any  effect  on  the  com¬ 
puted  velocities  profiles,  as  can  be  seen  from  Figure  11.  In  the  post-fracture-strength 
Equation  11,  the  strength  of  the  pulverized  ceramic  is  directly  proportional  to  the  applied 
pressure  when  a  =  0.  The  strength  is  zero  for  P  =  0;  twice  the  pressure  for  P  =  2,  and 
equal  to  the  pressure  for  p  =  1. 
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POST  FRACTURE  STRENGTH  EFFECT 


Figure  1 1 .  The  effect  of  post-fracture  strength  model  parameter 
P  on  the  computed  velocity  profOe. 


The  average  spall  signal  between  points  "E”  and  "F"  is  almost  the  same  for  these 
three  values.  It  appears  that  the  calibration  of  this  parameter  should  be  based  on  some 
other  type  of  impact  test  configurations  such  as  the  configuration  in  which  a  rod  pene¬ 
trates  into  a  ceramic  plate.  The  computed  stress  profiles  for  three  values  of  a  =  1,  5, 
and  10  kbars  did  not  influence  the  results  either;  therefore,  in  the  model  constant  calibra¬ 
tion  scheme,  a  was  set  to  zero  for  three  ceramics. 

Silicon  Carbide  (SiC) 

The  shock  portion  of  the  wave  profile  for  SiC  resembled  a  typical  wave  profile  of  a 
rate  sensitive  metal.  In  fact,  Kipp  and  Grady  (5)  successfully  modeled  the  low  and  high 
velocity  tests  on  SiC  using  a  strain  hardening  model  applicable  for  metals.  The  absence 
of  large  release  wave  dispersion  in  SiC  enabled  the  successful  use  of  a  simple  metal- 
based  model.  The  experimental  data  did  not  reveal  any  phase  transformation  in  SiC  up  to 
about  36  GPa.  The  HEL  kink  appears  at  about  15  GPa  in  both  tests. 

A  set  of  best  suitable  values  for  the  constants  were  determined  by  matching  the 
model  generated  velocity  profile  with  the  experimental  data  of  test  1  (see  Table  I).  The 
corresponding  values  are  given  in  Table  2.  The  equation-of-state  parameters  |3i,  ^2,  p3, 
and  r  are  given  in  Table  3.  These  values  were  employed  by  Johnson  and  Holmquist  (7) 
in  their  modeling  of  Kipp  and  Grady’s  experimental  data. 
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Table  2.  TTie  ceramic  model  constants 


Model  constants 

BaC 

SiC 

TiBz 

fo  (porosity  content) 

0.03 

0.0 

0.0 

Kic  (MPaVin)  (fracture  toughness) 

2 

4 

7 

p  (coefficient  of  friction) 

0.1 

0.1 

0.1 

N*  (m'^)  (number  of  flaws) 

5x10® 

1x10® 

1x10® 

"  1 

0.1 

0.1 

0.15 

ao  (microns)  (initiai  crack  size) 

0.5 

0.5 

2.0 

Table  3.  Equation-of-state  and  strength  constants 


Other  constants 

BaC 

SiC 

TiBa 

Density  (kg/m3) 

2516 

3177 

4452 

pi  (GPa) 

233 

204 

230 

Pz  (GPa) 

50 

204 

930 

p3(GPa) 

0 

0 

-1310 

r 

1 

1 

1 

Yhel  (GPa) 

12.5 

11.5 

12.0 

C3 

0.01 

0.01 

0.01 

G  (GPa) 

199 

187 

237 

Figure  12  compares  the  ceramic  model  generated  velocity  versus  time  plot  with  the 
Kipp  and  Grady  data  for  SiC  at  an  impact  velocity  of  1542  m/s.  The  model  reproduced 
all  the  salient  features  of  the  data.  As  the  next  step,  the  high-stress  experiment  (see  test 
2  in  Table  1)  at  an  impact  velocity  of  2100  m/s  was  simulated  using  the  same  set  of 
values  for  the  model  constants.  The  simulation  predicted  velocity  profile  matched  the 
data  extremely  well,  as  shown  in  Figure  13.  The  ceramic  model  reproduced  all  the 
salient  features  of  the  experimental  data. 

The  simulated  stress-strain  path  in  the  center  of  the  SiC  target  is  plotted  in  Figure  14. 
In  this  figure,  the  compressive  stress  and  compressive  pressure  are  plotted  as  positive 
quantities.  The  computed  HEL  is  about  18  GPa.  The  experimental  value  was  about 
15  GPa.  The  increased  value  in  the  computation  is  due  to  the  strain  rate  dependent  term^ 
in  Equation  7.  It  is  possible  to  lower  the  value  of  C|  and  adjust  the  value  of  C3  to 
match  the  computed  stress  at  HEL  with  the  data.  However,  for  the  sake  of  simplicity  in 
the  'model  constant  determination  scheme,  the  value  of  Y  that  corresponds  to  the  stress  at 
HEL  (see  Equation  8)  is  assigned  to  C|. 
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Rgure  12.  A  comparison  between  model  generated  velocity 
profile  with  the  Kipp  and  Grady  plate  impact  experimental  data 
at  impact  velocity  1542  m/s  for  silicon  carbide. 
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TIME  (MICROSECONDS) 


Figure  13.  A  comparison  between  model  generated  velocity 
profile  with  the  Kipp  and  Grady  plate  impact  experimental  data 
at  impact  velocity  2100  m/s  for  silicon  carbide. 
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Figure  14.  The  stress-strain  path  (under  one-dimensional  strain) 
in  the  midplane  of  an  SC  target  The  dotted  line  is  the  computed 
hydrostat 


Therefore,  the  transition  from  elastic  to  plastic  deformation  is  controlled  by  the  strain 
rate  dependent  strength.  A  value  of  12.5  GPa  was  employed  in  the  simulation.  Though 
Cl  was  chosen  to  be  equal  to  this  value,  the  high  e  (>IO'^/s)  slightly  increased  the  value  of 
Y.  This,  in  turn,  caused  a  slight  increase  in  the  computed  Kipp  and  Grady 

obtained  a  value  of  about  16  GPa  in  their  computational  analysis.  The  maximum  stress 
amplitude  in  compression  reached  about  26  GPa  and  elastically  unloaded  past  the  hydro¬ 
stat  (see  the  dotted  line  curve). 

At  about  0.02  strain,  the  stress  became  tensile  and  reached  a  tensile  maximum 
of  about  2  GPa.  The  ceramic  model  reduced  the  stress  amplitude  in  tension  due  to 
microcracking.  The  ability  to  control  the  tensile  stress  amplitude  is  an  important  feature 
of  the  model.  The  degradation  of  the  tensile  strength  due  to  microcracking  was  a 
continuous  process  (not  an  instantaneous  process)  and,  therefore,  the  stress  relaxation 
occurred  continually  until  the  microcrack  density  reached  a  critical  value  of  9/16,  accord¬ 
ing  to  the  Budiansky  and  O’Connell  solution  (19)  for  tensile  pressure  loading  conditions. 

Johnson  and  Holmquist  (7),  Steinberg  (6),  and  Addessio  and  Johnson  (8)  successfully 
reproduced  the  measured  velocity  profiles  of  Kipp  and  Grady  (5).  In  general,  the  shock 
behavior  of  SC  is  relatively  staightforward  compared  to  B4C  and  TiB2.  These  two  high 
strength  ceramics  exhibited  heterogeneous  deformations.  The  strain  rate  and  pressure 
hardening  behaviors  in  these  ceramics  did  not  resemble  the  metal  behavior. 
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The  pressure  versus  volumetric  strain  plot  showed  different  loading  and  unloading 
paths.  Since  the  model  degraded  the  bulk  modulus  according  to  microcracking  under 
tensile  loading,  the  unloading  path  (see  Figure  14)  between  the  strains  of  0.02  and  zero 
differed  from  the  loading  path.  The  slight  difference  between  the  loading  and  unloading 
paths  at  strains  above  0.04  was  due  to  microcracking  under  compression.  As  mentioned 
earlier,  the  crack  opening  is  possible  under  compressive  pressure  when  one  of  the  princi¬ 
pal  stresses  is  positive  (tension). 

The  time  histories  of  damage  and  axial  stress  at  the  midpoint  of  the  ceramic  target 
are  plotted  in  Figure  15.  To  compare  the  stress  with  damage,  the  stress  is  normalized 
by  the  Ohel  GPa)-  The  microcracking  occurred  at  point  "A"  under  compressive 

loading  and  reached  a  maximum  of  about  0.2  at  point  "B"  in  the  figure.  About  t  = 

1.12  ps,  the  stress  (see  the  dotted  line)  became  tensile  (negative  in  the  figure)  and  dam¬ 
age  reached  one.  In  the  model,  the  ceramic  material  at  the  midpoint  of  the  target  failed 
due  to  tensile  cracking  and  could  no  longer  support  any  tensile  pressure. 


OAmCE  AM)  STRESS  HISTORIES 
AT  TVC  MIDPOIMT  OF  TVC  TARGET 


Figure  15.  The  computed  damage  and  stress  histories  at 
a  midsection  of  the  SC  target.  The  axial  stress  is  normalized 
byoHEL=  15GPa. 


The  maximum  amplitude  of  the  tensile  stress  is  controlled  by  the  degradation  of 
shear  and  bulk  moduli  under  tensile  loading.  When  the  crack  density  reaches  a  critical 
value  of  9/16  in  the  Budiansky  and  O’Connell  (19)  relationship  (see  Appendix  A),  corre¬ 
spondingly  the  damage  reaches  one.  At  this  time,  both  the  moduli  become  zero.  In 
the  simulation,  the  element  has  failed  under  tension  and  can  no  longer  transfer  tensile 
pressure.  The  stresses  and  pressure  are  set  to  zero  for  the  elements  where  the  tensile 
damage  reached  the  value  of  one  at  point  "C".  The  elements  are  allowed  to  sustain 
compressive  stresses  until  the  element  pulverizes  under  compressive  crack  growth. 
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Boron  Carbide  (B4C) 


Figure  16  compares  the  computed  velocity  profile  with  the  experimental  measurement 
for  B4C  at  1546  m/s.  Kipp  and  Grady  (5)  and  Steinberg  (6)  provided  some  explanations 
for  the  unusual  behavior  of  this  material  under  high  velocity  impact  loading.  According 
to  Kipp  and  Grady,  the  sporadic  nature  of  the  wave  profile  at  the  HEL  is  attributed  in 
part  to  a  heterogeneous  failure  or  faulting  mechanism  during  the  initial  shock  loading. 

The  plateau  at  the  HEL  also  means  a  slower  shock  speed  compared  to  the  elastic  wave 
speed.  In  the  present  study,  an  attempt  is  made  to  model  the  overall  (or  average) 
response  of  B4C  and,  accordingly,  the  model  constants  were  adjusted  to  reproduce  the 
data  in  the  best  possible  manner.  The  best  constants  for  the  ceramic  model  to  describe 
the  complex  behavior  of  B4C  were  obtained  by  matching  the  computed  velocity  histories 
with  the  data  from  test  3  (see  Table  1).  The  corresponding  comparison  between  the 
model  and  the  data  is  given  in  Figure  16.  The  equation-of-state  parameters  and  the 
model  constants  for  B4C  were  given  earlier  in  Table  3  and  Table  2,  respectively. 

The  measured  velocity  history  in  the  high-stress  amplitude,  test  4  was  predicted  using  the 
same  set  of  constants. 


TIME  (MICROSECONDS) 


Figure  16.  Comparison  between  model  and  Kipp  and  Grady  (5) 
data  for  a  boron  carbide  target  (see  test  3  in  Table  1 ). 


Figure  17  compares  the  computed  velocity  history  with  the  measured  data  at  2210 
m/s.  The  model  matched  the  experimental  data  very  well.  The  model  could  not  cap¬ 
ture  the  entire  dispersive  nature  of  the  unloading  portion  of  the  velocity  profile;  how¬ 
ever,  the  damage  evolution  through  Uie  time  dependent  microcracking  process  led  to 
the  degradation  of  elastic  moduli  and  this  degradation  allowed  some  sort  of  dispersion 
of  the  unloading  wave.  In  general,  the  ceramic  model’s  ability  to  reproduce  the  data 
from  the  low  and  high  velocity  tests  of  B4C  is  quite  impressive. 
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Figure  17.  Comparison  between  model  and  Kipp  and  Grady  (5) 
data  for  a  boron  carbide  target  (see  test  3  in  Table  1). 


Titanium  Oiboride  (TlBa) 

This  material  typically  revealed  two  breaks  in  the  loading  wave  portion  of  the  meas¬ 
ured  profiles.  In  the  present  effort,  no  attempt  is  made  to  model  the  two  breaks.  The 
measured  velocity  histories  from  the  two  plate  impact  experiment;  one  at  1SS6  m/s  and 
the  other  at  207S  m/s,  are  compared  in  Figure  18.  The  initial  compressive  loading  in 
the  lower  velocity  experiment  (test  5)  is  very  dispersive  with  a  larger  rise  time  compared 
to  the  higher  velocity  experiment  (test  6).  The  release  wave  from  the  flyer  arrives 
sooner  in  the  higher  velocity  case.  Since  the  flyer  plate  thickness  in  test  6  is  smaller 
(3.337  mm)  compared  to  the  flyer  plate  in  test  5  (3.972  mm),  the  pulse  duration  is 
shorter  in  test  6.  In  addition,  the  computational  analysis  of  Kipp  and  Grady  indicates 
shock  wave  formation  above  13.7  GPa  in  TiBa.  This  will  also  lead  to  a  faster  shock 
wave  speed  in  the  high  velocity  test  compeu'ed  to  the  dispersive  loading  wave  in  the 

lower  velocity  test  and  resulted  in  a  shorter  pulse  duration. 

The  best  suitable  model  constants  were  determined  for  TiBa  by  matching  the  com¬ 
puted  velocity  history  with  the  data  from  test  5.  These  constants  were  given  in  Table  2. 
Figure  19  successfully  compares  the  computed  profile  with  the  data.  The  ceramic  model 
reproduced  all  the  salient  features  of  the  measured  velocity  profile  including  the  disper¬ 
sive  nature  of  the  loading  as  well  as  the  unloading  waves.  The  high  velocity  test  was 

predicted  using  the  set  of  constants  determined  by  matching  the  lower  velocity  test;  the 
comparison  is  shown  in  Figure  20.  The  model  showed  more  dispersion  of  the  loading 
wave  compared  to  the  experiment.  Assuming  that  TiB2  shocks-up  above  a  shock  stress 
of  about  32  GPa,  the  equation-of-state  in  conjunction  with  the  constitutive  equations  was 
unable  to  model  the  high  amplitude  experimental  data.  However,  a  decent  matching  can 
be  seen  between  the  spall  portions  of  the  computed  and  experimental  profiles. 
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TITANIUM  DIBORIDE  TEST  DATA 


TIME  (MICROSECONDS) 


Figure  18.  The  measured  velocity  profiles  from  the  two  plate 
experiments  of  Kipp  and  Grady  (5)  onTiBa  targets  are  compared. 


TIME  (MICROSECONDS)' 


Rgure  19.  A  comparison  between  model  generated  velocity  profile 
with  the  Kipp  and  Grady  plate  impact  experimental  data  at  impact 
velocity  of  1515  m/s  forTiBa. 
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Figure  20.  A  comparison  between  model  predicted  vetocttyproffle 
with  the  Kipp  and  Grady  plate  knpact  experimental  data  at  impact 
velocity  of  2100  m/s  for  TiB2. 


Summary  and  Conclusions 


Summary 

Recently,  Rajendran  (1,4)  reported  the  development  and  application  of  an  advanced  ce¬ 
ramic  model.  This  ceramic  model  has  been  implemented  into  the  1991  version  of  the  EPIC- 
2  code  by  Grove  (17)  and  successfully  used  to  model  the  impact  behavior  of  AD85  ceramic. 

There  are  eight  microcrack  model  constants:  N*,  a©,  p,  "2’  ”1’  "2 

A  porous  ceramic  requires  one  additional  constant,  fo,  which  is  the  initial  porosity.  The 
fracture  toughness  value  is  taken  from  fracture  mechanics  handbooks.  The  coefficient  of 
friction  p  is  assumed  to  be  equal  to  0.1  based  on  the  sensitivity  study.  This  value  is 
arbitrarily  assumed  for  SiC,  B4C,  and  TiB2.  The  tensile  crack  growth  factor,  n|,  and 
the  tensile  crack  growth  index,  n^,  are  assumed  to  be  equal  to  one;  this  eliminates 
the  requirement  for  determining  these  two  constants;  therefore,  the  ceramic  model  effec¬ 
tively  requires  only  four  microcrack  constants  and  one  porosity  constant. 

The  rest  of  the  constants  were  calibrated  by  matching  the  simulated  velocity  or  stress 
profiles  with  the  measure  profiles  in  a  plate  impact  experiment.  The  strength  model 
constants  Ci  is  assumed  to  be  equal  to  the  compressive  strength  value  at  HEL  (Yhel). 
The  rate  sensitivity  constants  C3  is  calibrated  by  matching  the  shape  of  the  elastic-plastic 
transition  portion  of  the  shock  wave  (velocity  or  stress)  profile. 

In  summary,  this  report  presented  the  impact  behavior  modeling  of  three  ceramics; 

SC,  B4C,  and  TiB2.  The  main  objective  of  the  present  work  was  to  determine  the  model 
constants  for  the  model  using  the  plate  impact  experimental  measurements  of  Kipp  and 
Grady  (5).  The  following  constants  are  assumed  to  be  the  same  for  all  three  ceramics: 
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(12) 


ji.  =  0. 1,  n|  =  1.0,  n  j  =  1.0,  and  n .,  =  1.0  . 

The  other  constants,  a©,  N*,  ni-,  f©,  and  C3  were  adjusted  to  give  the  best  agreement 
with  the  data  from  the  plate  impact  tests  at  1.5  km/s. 

Conclusions 

In  general,  the  inelastic  deformation  in  ceramic  materials  due  to  impact  loading, 
consists  of  elastic,  plastic,  and  microcracking  components.  It  is  not  experimentally  possi¬ 
ble  to  isolate  these  deformations  through  any  direct  measurements.  With  lack  of  any 
recovery  techniques  to  examine  the  post-impacted  ceramic  targets,  especially  at  high  veloc¬ 
ity  impact,  speculative  assumptions  can  only  be  made  on  the  various  deformation  and 
fracture  processes. 

The  computational  analysis  of  the  impact  experiments  using  advanced  ceramic  models 
will  indeed  help  in  evaluating  the  various  possible  deformation  and  fracture  modes  in  the 
ceramic  materials.  So  far  the  experimentalists  have  been  making  accurate  velocity  and 
stress  measurements,  high  speed  photographs,  and  X-ray  radiographs.  Interpretations  and 
validations  of  these  valuable  measurements  demand  a  detailed  computational/analytical 
modeling  of  the  impact  experiments.  These  modeling  efforts  will  eventually  lead  to  a 
greater  understanding  of  the  impact  behavior  of  ceramics. 

In  order  to  increase  confidence  in  using  the  ceramic  models  as  predictive  tools  in 
armor/anti-armor  applications,  generality  of  the  model  parameters  should  be  tested  using  a 
variety  of  experimental  data.  The  experiments  should  not  only  include  various  stress- 
strain  states  but  also  a  range  of  velocity  regimes.  Though  the  high  velocity  (2.2  km/s) 
plate  impact  experiment  could  well  be  successfully  predicted  using  the  constants  deter¬ 
mined  from  the  low  velocity  (2.5  km/s)  one-dimensional  experiment,  the  generality  of  the 
model  constants  required  further  validation. 

For  this  purpose,  other  impact  test  configurations  with  different  stress  states  should 
be  modeled  using  the  preliminary  constants  presented  in  this  report.  However,  availabil¬ 
ity  of  a  preliminary  set  of  model  constants  for  armor  ceramics  will  encourage  the  scien¬ 
tists  and  engineers  who  perform  advanced  computational  calculations  in  their  design 
studies  of  armor/anti-arnor  to  use  advanced  material  constitutive  and  failure  models. 

In  conclusion,  this  microphysical  ceramic  model  is  capable  of  modeling  the  impact 
behavior  of  the  intact  ceramic  up  to  pulverization.  Once  the  ceramic  is  pulverized,  the 
accurate  and  realistic  description  of  the  confined  pulverized  ceramic  requires  improved 
equation-of-state  and  impact  strength  behavior  models.  In  the  computer  code  calculations, 
a  simplistic  Mohr-Coulomb  type  strength  model  (see  Equation  II)  is  usually  employed  to 
describe  the  powdered  ceramic.  While  the  microphysical  ceramic  model  is  going  to  cer¬ 
tainly  improve  the  quality  of  computational  armor/anti-armor  design  calculations,  the 
absence  of  a  microphysical  model  for  the  fractured  (pulverized)  ceramic  unfortunately 
diminishes  the  ceramic  model  capability  in  predicting  the  depth  of  penetration/perforation 
of  projectiles.  Presently,  Namat-Nasser  and  his  students  at  the  University  of  California  in 
San  Diego,  CA,  Clifton  and  his  students  at  Brown  University  in  Providence,  RI,  Curran 
and  his  coworkers  at  SRI  International,  and  several  other  researchers  are  making  progress 
towards  the  development  of  dilatational  constitutive  models  for  powdered  ceramics. 
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Appendix 


Degraded  Moduli  Expressions 

The  elastic  stress-strain  relationship  in  Equation  2  corresponds  to  the  microcracked 
ceramic  material.  The  elements  of  the  stiffness  matrix  M  in  Equation  2  is  obtained  from 
the  derivations  of  Margolin  (11).  Margolin  derived  the  degraded  elastic  moduli  for  nonin¬ 
teracting,  penny-shaped  microcracks  of  various  sizes  and  in  random  orientations.  The  com¬ 
ponents  of  the  compliance  (moduli)  tensor  Cijici  which  is  the  inverse  of  the  stiffness 
tensor  Mijkl  for  the  isotropic  elastic  moduli  in  Equation  3  is  given  by; 


Qjkl  =  Cl  6ik5ji  +  C2  5ii5jk  +  CsSjjSki 
where 


(13) 


1 

Cl  -  Bq  +  , 

C2  =  Dq  +  . 


and 


(14) 

(15) 


C3  ~  Aq 

(16) 

2( 1 +v)G 

In  the  above  equations,  G  and  v  are  the  shear  modulus  and  Poisson’s  ratio,  respectively, 
of  the  undamaged  material  while  Aq,  Bq,  and  Dq  are  damage  parameters  whose  values 
depend  on  the  stress  state.  To  evaluate  these  parameters,  Margolin  defined  a  microcrack 
density  parameter  as 


16N: 


a 


3 

max 


45  E 


(17) 


In  Equation  17,  y*  is  the  microcrack  density  and  E  is  the  Young’s  modulus  of 
the  undamaged  material.  N*  is  the  number  of  microcracks  per  unit  volume,  and  amax 
is  the  maximum  microcrack  size.  Margolin  identified  the  following  four  cases  of  stress 
state  in  evaluating  the  damage  parameters  Aq,  Bo,  and  Do'. 


Case  j: 


(ail  prinicpal  stresses  are  tensile) 


Ao  =  (1  -v^)-(l  +v)  Y* 


(18) 
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Bo  =  (1  -v-)  +  4(l  +v)  y 


(19) 


Dq  —  Aq 


(20) 


Case  2: 

o, ,02.03  <  0 

(all  principal  stresses  are  compressive) 


Ao  =  -(1  +v)Y* 


(21) 


Bo  =  4(1  +v)Y* 

(22) 

Dq  —  Aq 

(23) 

Equations  21  through  23  will  result  in  the  degradation  of  the  shear  modulus  but  not  of 
the  bulk  modulus  because,  under  compression,  only  crack  movement  of  the  closed 
microcracks  under  modes  II  and  III  are  permitted. 

Case  3: 

,  o’2 ,  >  0 ,  ^3  <  0 

(two  principal  stesses  are  tensile  and  one  principal  stress  is  compressive) 


Aq  — 


( 5B^  -  36^)  2 

^  (1-v^)  -  (1+v) 

2 


Bo  = 


(1  -v^)  +  4(1  +v) 

2 


Do  =  [(6p5-5p3)  (l-v2)  -  (l+v)]  y* 


where 


P  = 


o,  +  03-203 


(24) 

(25) 


(26) 


(27) 
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Case  4: 


o,  >  0.  a^.Oj  <  0 

(one  principal  stress  is  tensile  and  two  principal  stresses  are  compressive) 


Ao  = 


Bo  = 


5(1 -P^)- 3(1 -8^) 


(1  -  V  "■)  -t-  4  (1  +v  ) 


Do  =  I[6(l-p5)-5(l-|3^)](l-v2)  -  (1+v)  1  Y* 


(28) 


(29) 


(30) 


where 


- 


°2  °3 


-  2<J  j 


(31) 


In  the  ceramic  model,  the  compliance  tensor  C  is  analytically  inverted  to  the  stiffness 
tensor  M  using  the  following  identity  relationship: 


Mr-  Sji  +  8ii  5jk ) 
^  ijmn  ^  mnkl  2 


(32) 


A  relationship  can  be  established  between  Cj,  C2,  and  C3  of  Equations  14  through  16 
and  mi,  m2,  and  m3  of  the  following  relationship: 


M  jjij,  =  mi  5ik  5ji  +  m2  5ii  6jk  +  m3  Sjj  5ki 

where 

1 


(33) 


mi  =  m2  = 


and 


m3  = 


2  (C,  +  C2) 


-C3 


(34) 


(Cl  +  C2)(C]  +  C2  +  3C3) 


(35) 
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The  effective  degraded  shear  and  bulk  moduli  are  defined  as 


G  =  m 


and 


(36) 


K 


=  nt3  +  2  mi 


(37) 


Instead  of  employing  Margolin’s  expression  for  Case  1  (all  principal  stresses  are  ten¬ 
sile),  Budiansky  and  O’Connell  (19)  solutions  can  be  used  for  randoin!,  oriented,  non¬ 
interacting  microcracks  under  tensile  loading.  The  corresponding  relationships  for  the 
damaged  stiffness  solutions  are: 


mi  =  m2  =  G 


(38) 


m3 


2G  V 
(1-2V) 


(39) 


where 


( 

V  =  V 

s 

and 


(40) 


G  =  G 


f'32'l(r  -  v)(5  -  V)  1 
[45j  a-v) 


(41) 


In  these  equations,  v  and  G  are  the  Poisson’s  ratio  and  shear  modulus,  respec¬ 
tively,  of  the  microcrack  damaged  material.  Using  m|  and  m3  from  Equations  38  and 
39,  the  degraded  bulk  modulus  can  be  computed  from  Equation  37.  It  is  obvious  from 
Equations  38  through  41  that  a  complete  loss  of  strength  is  predicted  when  the 
microcrack  density  y  reaches  9/16.  For  tensile  loading  conditions,  based  on  the  compari¬ 
son  with  Margolin’s  equations,  there  is  no  bound  on  the  crack  density;  however, 
Budiansky  and  O’Connell’s  solution  limits  the  crack  density  to  9/16.  This  permits  the 
damage  parameter  to  vary  from  zero  (no  damage)  to  one  (fully  damaged);  therefore,  in 
the  ceramic  model,  Budiansky  and  O’Connell’s  equations  are  used  instead  of  Margolin’s 
equations  for  the  case  when  all  the  principal  stresses  are  positive  (see  Case  1). 
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